Solving for the dynamics of the universe 

Valerio Faraoni 

Research Group in General Relativity (RggR) 
Universite Libre de Bruxelles 
Campus Plaine CP 231, Blvd. du Triomphe 
1050 Bruxelles, Belgium 
E-mail: vfaraoni@ulb.ac.be 



Abstract 

A new method of solving the Einstein-Friedmann dynamical equations of a spa- 
tially homogeneous and isotropic universe is presented. The method is applicable 
when the equation of state of the material content assumes the form P = (7 — 1) p, 
7 = constant. The solution for the expansion factor is commonly given only for 
7 = 0, 1 and 4/3 when the curvature index is K = ±1. The proposed procedure 
is valid for general 7 and K and it allows for ease of derivation of the solutions. 
This alternative method is useful pedagogically to introduce basic cosmology. 
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1 Introduction 



The Friedmann-Lemaitre-Robertson-Walker (FLRW) solutions to the Einstein field 
equations of general relativity are a cornerstone in the development of modern cos- 
mology. The FLRW metric describes a spatially homogeneous and isotropic universe 
satisfying the Copernican principle, and is the starting point for studying dynamical 
models of the universe. To this end, one must specify the nature of the matter which is 
the source of gravitation by assigning the corresponding equation of state, which varies 
during different epochs of the history of the universe. 

The FLRW metric is given, in comoving coordinates (t, r, 6, tp), by the line element 

ds 2 = -dt 2 + a 2 (t) [dr 2 + f 2 (r) (d6 2 + sin 2 6dip 2 )] , (1.1) 

where the three possible elementary topologies are classified according to a normalized 
curvature index K which assumes the values 0, ±1, and 

/(r) = sinhr (K = -1) (1.2) 

for the open universe, 

f(r) = r (K = 0) (1.3) 

for the critical universe, and 

/(r) = sinr (K = +1) (1.4) 

for the closed universe. Other topologies are possible (see e.g. Ref. 1, p. 725). The 
function a(t) of the comoving time t (the "scale factor") is determined by the Einstein- 
Friedmann dynamical equations 2-7 

a AnG , , 
- = - — (p + 3P), (1.5) 

-) 2 ^-^, (1.6) 

where p and P are, respectively, the energy density and pressure of the material content 
of the universe, which is assumed to be a perfect fluid. An overdot denotes differentiation 
with respect to the comoving time t, G is Newton's constant and units are used in which 
the speed of light in vacuum assumes the value unity. It is further assumed that the 
cosmological constant vanishes. 
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One can solve the Einstein-Friedmann equations ( |1 . 5| ) and ( |1.6| ) once a barotropic 
equation of state P = P(p) is given. In many situations of physical interest, the equation 
of state assumes the form 

P = (7 — 1) p , 7 = constant . (1.7) 

This assumption reproduces important epochs in the history of the universe. For 7 = 1 
one obtains the "dust" equation of state P = of the matter-dominated epoch; for 
7 = 4/3, the radiation equation of state P = p/3 of the radiation-dominated era; for 
7 = 0, the vacuum equation of state P = — p of inflation. For 7 = 2/3, the curvature- 
dominated coasting universe with a(t) = a^t is obtained. 

When the universe is spatially flat (K = 0), the integration of the dynamical equa- 
tions (|1.5|) , (|1.6| ) with the assumption (|1.7|) is straightforward 2 for any value of the 
constant 7: 

a{t)=a t^ (7^0), (1.8) 

a (t) = a e m , H = (7 = 0). (1.9) 

However, for K — ±1, the solution is given (explicitly or in parametric form) only for 
the special values 7 = 1,4/3 in old and recent textbooks 1-7 . 

It is actually not difficult to derive a solution for any nonzero value of the constant 
7 using a standard procedure 3 which is general, i.e. it does not depend upon the as- 
sumption P = (7 — 1) p. Section 2 presents a summary of the usual method for deriving 
the scale factor a(t). An alternative method, which consists in reducing the Einstein- 
Friedmann equations to a Riccati equation, and is applicable when Eq. ( |1.7|) holds, is 
explained in Sec. 3. From the mathematical point of view, the latter method avoids the 
consideration of the energy conservation equation 

p + 3(P + p)- = (1.10) 
a 

and the calculation of an indefinite integral; rather, one solves a nonlinear Riccati equa- 
tion. The alternative approach is more direct than the standard one and the gain 
in clarity of exposition makes it more suitable for an introductory cosmology course. 
Section 4 presents a brief discussion of the equation of state P = (7 — 1) p and the 
conclusions. 
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2 The standard derivation of the scale factor 



The standard method to obtain the scale factor a of the universe proceeds as follows 3 . 
The energy conservation equation (|1.10|) yields 

[dp 

31na = — / h constant (2.1) 

J P + p 

for 7 7^ 0. Upon use of the conformal time 7] defined by 

dt = a(r])dr] , (2.2) 
the Einstein-Friedmann equations (|1.5| ), ( [Of ) yield 



da (2 3) 



In order to obtain a(rj), one prescribes the equation of state P = P(p), solves Eq. fl2.1|) 
and inverts it obtaining p = p(a). Further substitution into Eq. ( |2.3j ) and inversion 
provide a = a(rj). Integration of Eq. (|2.2j ) provides the comoving time t(rj) as a function 
of conformal time. The scale factor is then expressed in parametric form (a(j]),t(j])). 
Sometimes it is possible to eliminate the parametric dependence on r\ and obtain the 
expansion factor as an explicit function a(t) of comoving time. 

The method is quite general; as a particular case, it can be applied when the equation 
of state is of the form P = (7 — 1) p with constant nonvanishing 7. Equations (|2.1| ) and 
(2.3) then yield 

a^p lh = constant , (2.4) 
y = ±f , da (2.5) 



3 

for 7 7^ 0, 2/3, where C\ is an integration constant. By introducing the variable 

x= (*^lf*a, (2.6) 



and using 



/ 



= i l JV^ z A (27) 

xyjx n + 1 n \ Vx n + 1 + 1 J 
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dx 



one integrates and inverts Eq. 



x\/x n — 1 
TBT) to obtain 



2 

— arcsec 
n 



(x n ' 2 ) 



a(rj) = a sinh 1//c (c?7) 



for K = — 1. a and c are constants, with 



rv 

/ <ir/ sinh icq') 
Jo 



7-1 



(2.8) 

(2.9) 
(2.10) 

(2.11) 



and the boundary condition a (77 = 0) = has been imposed. 
Similarly, for K = +1, one obtains 



a(r/) = a [cos (cr? + d)] 1 ^ 
rv 

t(rj) — Oo / g??7 [cos (077' + d)] 



(2.12) 
(2.13) 



l/c 



For 7 = 2/3 and K — — 1 one obtains a curvature-dominated universe for which Eq. ( |1.6|) 
is approximated by (a/a) 2 ~ —K/a 2 . In this case, Eq. fl2.3| ) yields a = aoexp(/3?7) 
(/3 = constant), and t = toe'" 3 gives a = a$t. It is easier to obtain this form of the scale 
factor directly from Eq. ( |1 . 5[ ) , which reduces to a = 0. 

The solutions ( |2.9| ), ( |2.10| ) and ( |2.12| ), ( |2.13|) for the scale factor are presented in the 
textbooks only for the special values 1 and 4/3 of the constant 7. For 7 = 4/3 one 
eliminates the parameter 7] to obtain 



for K = +1, and 



a(t) = a 



a(t) = ao 



1 + 



21 1/2 



1/2 



(2.14) 



(2.15) 



for K = -1. 

The standard solution method of the Einstein-Friedmann equations has the virtue of 
being general; it does not rely upon the assumption (|1.7f) . However, the need to invert 
Eqs. ( |2.1| ) and (|2.3|) and to compute the indefinite integrals ( |2.6| ), ( 2.7|) detracts from the 
elegance and clarity that is possible when the ratio P/p is constant. The latter condition 
is satisfied in many physically important situations. 
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3 An alternative method 

There is an alternative procedure to derive the scale factor for a general value of 7 when 
the equation of state of the universe's material content is given by P = (7 — 1) p with 
7 = constant, which covers many cases of physical interest. Being straightforward, this 
new method is valuable for pedagogical purposes and proceeds as follows: Eqs. ( |1.5|) , 
(|OJD and (Q yield 

a ( a \ 2 cK 

~ + c - + — = > 3- 1 

a \aj er 

with c given by Eq. fl2.11| ). For K — 0, Eq. (|3.1| ) is immediately integrated to give 
Eqs. flDj), QO|) . For K = ±1, Eq. Q is rewritten as 

— + (c-l) f-") +c^ = 0, (3.2) 
a \ a / 

by making use of the conformal time 77, and where a prime denotes differentiation with 
respect to r\. By employing the variable 

a' 

u = — , 3.3 

a 

Eq. ( |3.2|) becomes 

u' + cu 2 + Kc = , (3.4) 
which is a Riccati equation. The Riccati equation, which has the general form 

dy 

— = a(x)y 2 + b(x)y + c(x) (3.5) 
ax 

where y = y(x), has been the subject of many studies in the theory of ordinary differential 
equations, and can be solved explicitly 8,9 . The solution is found by introducing the 
variable w defined by 

«=i^, (3.6) 
c w 

which changes Eq. ( |3.4|) to 

w" + Kc 2 w = , (3.7) 

the solution of which is trivial. For K — +1, one finds the solutions Q2.12| ), ( [2.13| ), while 
for K = — 1, one recovers (|2.9| ), ( j2.10| ). 
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The alternative method is applicable when the equation of state assumes the form 
P = (7 — 1) p and the cosmological constant vanishes. With these conditions satisfied, 
the alternative solution procedure is more direct than the general method. It is now 
appropriate to comment on the equation of state P = (7 — 1) p, 7 = constant, which 
reproduces several situations of significant physical interest. 



4 Discussion 

The assumption that the equation of state is of the form P — (7 — 1) p with constant 7 is 
justified in many important situations which describe the standard big-bang cosmology. 
However, it is important to realize that Eq. ( |1.7| ) is a strong assumption and by no means 
yields the most general solution for the scale factor of a FLRW universe. To make a 
physically interesting example, consider the inflationary (i.e. a > 0) epoch of the early 
universe in the K = case. Many inflationary scenarios are known 10 , corresponding to 
different concave shapes of the scale factor a(t); the assumption (|1.7|) allows the solutions 

a(t) = a t^ (4.1) 

("power-law inflation") for < 7 < 2/3, and 

a = a e m , H = (4.2) 

for 7 = 0. Vice-versa, using the dynamical equations ( |1.5|) , (|1.6| ), it is straightforward to 
prove that the latter solutions imply P/p = constant. The assumption ( |1.7| ) reproduces 
only exponential expansion (the prototype of inflation) and power-law inflation. All 
the other inflationary scenarios correspond to a j(t) which changes with time during 
inflation. A time-dependent 7(i) can also be used to describe a non-interacting mix- 
ture of dust and radiation; however the method of solution of the Einstein-Friedmann 
equations ( |1.5| ), ( |1.6| ) presented in Sec. 3 applies only when 7 is constant, and when the 
cosmological constant vanishes. If the cosmological constant is nonzero, Eq. (|3~T| ) does 
not reduce to a Riccati equation. 

The limitations of the assumption 7 = constant are thus made clear: while this new 
approach to the Einstein-Friedmann equations does not replace the standard approach, 
it is more direct and is preferable in an introduction to cosmology. Its value lies in the 
ease of demonstration of the solutions, which is crucial for students to grasp the basic 
concepts of cosmology. 
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